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Abstract 

The light-front (LF) canonical quantization of quantum chromodynamics in covariant 



> 

m 

gauge is discussed. The Dirac procedure is used to eliminate the constraints in the gauge- 
Q\ ■ fixed front form theory quantum action and to construct the LF Hamiltonian formulation. 

The physical degrees of freedom emerge naturally. The propagator of the dynamical tp + part 
of the free fermionic propagator in the LF quantized field theory is shown to be causal and 
not to contain instantaneous terms. Since the relevant propagators in the covariant gauge 
formulation are causal, rotational invariance — including the Coulomb potential in the static 
limit — can be recovered, avoiding the difficulties encountered in light-cone gauge. The Wick 
rotation may also be performed allowing the conversion of momentum space integrals into 
Euclidean space forms. Some explicit computations are done in quantum electrodynamics to 
illustrate the equivalence of front form theory with the conventional covariant formulation. 
LF quantization thus provides a consistent formulation of gauge theory, despite the fact that 
the hyperplanes x^ = used to impose boundary conditions constitute characteristic surfaces 
of a hyperbolic partial differential equation. 
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1 Introduction 

The quantization of relativistic field theory at fixed light-front time r = (t — z/c)/\^2, 
which was proposed by Dirac [I] half a century ago, has found important applications 
DH, H, £| in gauge theory and string theory ||. The light-front (LF) quantization of QCD 
in its Hamiltonian form provides an alternative approach to lattice gauge theory for the 
computation of nonperturbative quantities, such as the spectrum and the light-cone 
Fock state wavefunctions of relativistic bound states. LF variables have found natural 
applications in several context, for example, in the quantization of (super-) string 
theory and M-theory ||. It has also been employed in the nonabelian bosonization || 
of the field theory of N free Majorana fermions. 

Since LF coordinates are not related to the conventional coordinates by a finite 
Lorentz transformation, the descriptions of the same physical result may be different 
in the equal-time (instant form) and equal LF-time (front form) of treatments of the 
theory This was also found to be the case in the recent study of some soluble two- 
dimensional gauge theory models, where it was also demonstrated that LF quantization 
is very economical in displaying the relevant degrees of freedom, leading directly to a 
physical Hilbert space. LF-time-ordered perturbation theory has also been applied || |9j 
to massive fields. It was used in the analysis of the evolution of deep inelastic structure 
functions [ITU] , the evolution of the distribution amplitudes which control hard exclusive 



processes in structure functions jL0|, and the evolution of the distribution amplitudes 
which control hard exclusive processes in QCD |Tl| . LF-time-ordered perturbation 
theory is much more economical than equal-time-ordered perturbation theory, since 
only graphs with particles with positive LF momenta p + = (p° +p 3 )/^/2 appear. It 
has been conventional to apply light-front quantization to gauge theory in light cone 
gauge A + = (A°+A 3 )/y/2 = 0, since the transverse degrees of freedom of the gauge field 
can be immediately identified as the dynamical degrees of freedom, and the ghost fields 



can be ignored in the quantum action of the nonabelian gauge theory ||12| . However 



it does bring severe penalties, such as the breaking of manifest rotational invariance, 
cumbersome LF time instantaneous interactions, an analytically complex gauge-field 
propagator, and a difficult renormalization procedure [DJ . In addition, light-cone gauge 



complicates obtaining the Coulomb limit when one applies the LF formalism to simple 
nonrelativistic systems. 

In this paper we will discuss the LF quantization of (massless) gauge field theory 
in covariant gauges such as Feynman gauge. The penalty in reintroducing Faddeev- 
Popov ghosts and other fields with ghost-like metric is compensated by the restoration 
of the Lorentz symmetries. Our general procedure will be first to study the gauge-fixed 



quantum action of the theory on the LF. The self-consistent Hamiltonian formulation 



is constructed following the Dirac method [JL4| for constrained systems. It is worth 
remarking that even after the addition of gauge-fixing and ghost terms in the front 
form theory, several second class constraints in general remain in the theory. They 
can be eliminated by constructing the Dirac brackets, and the theory can be quantized 
canonically by the correspondence principle in terms of a reduced number of independent 
fields. The commutation relations among them are also found by the Dirac method, 
and they are useful in order to obtain the momentum space expansions of the fields. 
The nondynamical projection of the fermion field will be eliminated using a nonlocal 
constraint equation. The Dyson- Wick perturbation theory expansion |TSJ based on LF 



coordinates then becomes straightforward Pj. In fact, the Feynman gauge has been 



already used [16[] in QED in the framework of Weinberg's time-ordered perturbation 
theory in the infinite-momentum-frame, where it was used to renormalize and compute 
the order a, a 2 and part of the a 3 radiative corrections to the electron anomalous 
moment. (In this gauge the theory on the LF resembles that of an interacting theory 
of four scalar fields {A + , A- 1 -, A~) with the matter fields. The covariant gauge thus 
should be feasible within the LF framework.) The challenge of renormalization in 



light-cone gauge and some of its aspects in Feynman gauge were addressed in |fP7 |. 
The front form theory has recently been used [18[ to quantize a nuclear physics model 
in which massive vector mesons couple to field-theoretic nucleons. In this approach, a 
canonical transformation || [| of the dynamical nucleon ip + field is used to simplify the 
constraint equation for ip_ (in abelian theory). The inclusion of fields with ghost-like 



metric in discretized light-cone quantization (DLCQ) |19| has been discussed in f20|j . 

The propagator for the dynamical component if>+(x) of the free fermionic field on 
the LF is shown to be causal, and it is demonstrated that, in the context of LF-time- 



ordered Dyson- Wick perturbation theory, it has no instantaneous term. This effect |TT 



is replaced now by explicit seagull and other interaction terms which can be evaluated 
systematically. Such interactions can be incorporated into a nonperturbative approach 
such as DLCQ. Next QCD is canonically quantized on the LF in covariant Feynman (or 
Landau) gauge. A novel feature of the gauge theory interactions in covariant gauges on 
the LF is the off-diagonal couplings of the A + and A~ fields. Some illustrations are given 
on the application the Dyson- Wick perturbation theory expansion based on the Wick 



theorem, built on the LF-time-ordered product. We recall that it was also used [21] to 



renormalize two-dimensional scalar field theory on the LF, with nonlocal interaction. 



It was also shown there |2T[ that even on the LF one may perform the Wick rotation 
in momentum space integrals and go over to the Euclidean-space integrals, rendering 
the computation of high order corrections as straightforward as in conventional theory. 



2 Notation 

The Lagrangian density corresponding to the quantum action \T1\ of QCD is described 
in standard notation by 

Cqcd = -\F a ^F\ v +B a d^A^ + ^B a B a +td^Xi a V a %X2 C +^ l ^D l ^-m^W (1) 
Here ipi is the quark field with color index j = l..N c for an SU(N C ) color group, A a 



/< 



the gluon field, Fj^ = d^A a v - d v A a ^ + gf abc A\A c u the field strength, V a % = (5 ac d^ + 
gf abc A\), D ij ^ = {8 ij d^-igA a ^(X a /2Y^, a = l..(N c 2 - 1) the gauge group index 
etc. The covariant gauge-fixing is introduced by adding to the Lagrangian the linear 
gauge- fixing term (B a d^A afl +^B a B a ) where B a is the Nakanishi-Lautrup auxiliary field 
and £ is a parameter. The xi a an d X2 a are the (hermitian) anticommuting Faddeev- 
Popov ghost fields, and the action is invariant under the BRS J23J transformation. 

The LF coordinates /i = +,—,1,2 are defined as x^ = (x + = X- = (x° + 
x 3 )/^/2, x~ = x + = (x° — x 3 )/\^2, x- 1 ) where x 1 - = (x 1 ^ 2 ) = (— Xi,— x%) are the 
transverse coordinates and \x = +, — , 1,2. The x + = r will be taken as the LF-time 
coordinate while x~ is the longitudinal spatial one. We can of course make the con- 
vention with the role of x + and x~ interchanged. The equal-x + quantized theory does 
seem to already contain the information on the equal-x~ commutator 0. For the four- 
momentum components we have k^ = (k~ = (k° — k 3 )/^/2, k + = (k° + k 3 )/\/2, Ar 1 ) 
where k + indicates the longitudinal momentum and k~ the LF energy. The LF compo- 
nents of any tensor, for example, the gauge field (A + , A~ , A^) are similarly defined, and 
the metric tensor g^ v may be read from A-B = A + B~ + A~ B + — A x B x . The 7 defined 
by 7 ± = (7 ± 7 3 )/v/2 satisfy (7+) 2 = (7") 2 = 0, 7%+ = 7-7 , 7+7" + 7-7+ = 21, 
7 + = 7~ and it follows that A^ = ^7 =F 7 ± = "757°7 ± are hermitian projection op- 
erators: (A ± ) 2 = A ± , A + A~ = A~A+ = 0, 7°A+ = A - 7°. The corresponding 
± projections of the LF Dirac spinor are ip± = A ±r ip and ip = ip + + ip-,ip = ip^ = 



3 Spinor field propagator on the LF 

The quark field term in LF coordinates reads 

+ $_ [zV^W.^J + (i^D ij ± - m5 ij )ip + j ] (2) 



where D l i± = (5^d± —igA a ±(X a /2)^). It shows that the minus components ip_ 3 are in 
fact nondynamical ( Lagrange multiplier ) fields without kinetic terms. The variation of 
the action with respect to ipi _ and ip J _ leads to the following gauge covariant constraint 
equation 

iV2D ij _i/;J = -(z 7 V^i - m-y°5 ij )i> + j , (3) 

and its conjugate. The ■0- 7 '_ components may thus be eliminated in favor of the dy- 
namical components ?/>+ 



ip- 3 (x) 



V2 



U-\x\A_)—U(x\A_) 



ijk 



(t^ ± D kl ± - m^°8 kl )iP + l {x). (4) 



Here, for a fixed r, U = U(x\A_) is an iV c x N c gauge matrix in the fundamental 
representation of SU(N C ) and it satisfies 



a_J7(s|A_) = -igU(x\A_) A_(z) 
with A_ = A a _X a /2. It has the formal solution 

U(x-,x^) = U(x- ,x^)Vex P {- W f X ~ d y-A^y-, X -) 



(5) 



(6) 



x~ 



where V indicates the anti-path-ordering along the longitudinal direction x~ . U has a 
series expansion in the powers of the coupling constant. 

The free field propagator for ip + is determined from the following quadratic terms 
in the Lagrangian (suppressing the color index) density 



iV2i/)+ t d+iJ) + + ^ + t (i7°7 ± 9 ± - m 7 °)^_ 



(7) 



contained in (2). Here we have used the free field constraint equation 2id_if}_ = 
(i r y L d± + m)^ + ip + which determines the dependent field ?/>_. The equation of motion 
for the independent component ip + is nonlocal in the longitudinal direction 



Ad + + (m + 27" L (9_l)7 Tr( m + i>l ®±)l 



Vm 



o. 



(8) 



The free field Hamiltonian formulation can be constructed by following the Dirac proce- 
dure []14j . The constraint equation arises now as a second class constraint on the canon- 
ical phase space. The Dirac bracket which takes care of these constraints is easily con- 
structed. The effective free LF Hamiltonian is found to be 7i l f = —ip + (i'~y ± d± —m)ip- 
and the canonical quantization performed by the correspondence of the Dirac brackets 



with the (anti-) commutators leads to the following nonvanishing local anticommutation 
relation 



{^+(r, x~, x 1 ), i/>+\t, y'^ 1 )} = ~7^ A+S ( X ' - y~)P{x L - y 1 ) 



(9) 



on the LF. They were proposed earlier in Ref. [8]. The equation of motion (8) for i\)+ 
is recovered as an Heisenberg equation of motion if we employ (9). 

The propagator in momentum space may be derived by going over to the Fourier 
transform of ip(x) over the complete set of linearly independent plane wave solutions 
of the free Dirac equation, say, for p + > 0. Such a set is spanned by uS r ' (p)e~ ip ' x 
together with v^ r '(p)e ip ' x where u^ r '(p) and v^ r '(p) = vS r >(p) c = C7 vS r '(p)* are linearly 
independent solutions of the Dirac equation in momentum space: (^ tJ, p fl —m)u <<r \p) = 
and p ■ x = (p~x + + p + x~ — p^ar 1 ). 

A useful form [33 of the solution for the four-spinors in the context of the LF 
quantization is 



u {r) (p) 



(vV 



171)2 



y/2p + A + + (m + <y ± p±)A~ 



fl W 



(10) 



ru^ with S 3 = «7 1 7 2 



where the constant spinors u^ satisfy j°u^ = u^ and T, 3 u^ 

and r = ±. The normalization and the completeness relations are: vf J '\p)u^ s \p) = 

S rs = -v {r \p)v {s \p), E r =±u {r) (p)u {r) (p) = (ft + m)/2m, E r =±v {r) (p)v (r) (p) = 

(ft — m)/2m and C is the charge conjugation matrix fl5|| . 

1 

Its A + projection is by construction very simple, vS r '+(p) = (\/2p + /m) 2 (A + u^) 
and they are eigenstates of S3 as well while the vS r ' correspond to the rest frame spinors 
for which \plp ± = m. 

The Fourier transform expansion may be written as 

1 



ip(x) 



(27T) 



r=± 



J2 / d 2 p ± dp + 6(p A 




b {r) (p)u {r) (p)e- ip - x + d Kr) (p)v (r) (p)e ip - x 



(11) 
The presence of the factor 6(p + ) in (11) follows also from the considerations of the 

covariant phase space (or LIPS) factor which is found relevant [2Jj| in the context of 

the analysis of the physical processes. For ip + = A + ip we find 



i[) + (x) 



(V2)i 



(2*0 



J2 f d 2 p ± dp + e( P + ) b ir \p)u i ; ) e- ip - x + d^ r \p) 



r=±- 



$(r) e ip.x 



(12) 



where u and v are constant spinors, and the integrations are to be taken from —00 
to 00 not only for p L but also over p + . It can be verified that the anticommutation 



relations (9) are satisfied if the creation and the annihilation operators are assumed to 
satisfy the canonical anticommutation relations, with the nonvanishing ones given by: 
{b^(p),tf {s \p')} = {d^(p),d^ s \p')} = 5 rs 8(p+ -p' + )5 2 (p^ -p' 1 ). 

The free propagator for the independent component ip+ is easily derived using (12) 
and (10) and it is given by 

(0\T(ij> +A (x)ij>l B (0))\0) = (0| [9(r)i; +A (x)iji B (0) - 9(-r)^ +B (0)^ +A (x)} |o) 

= ^^% J AW% + ) [9(T)e** - 9(-r)e^}(13) 

where A,B= 1.4 label the spinor components and we have used (12). The only relevant 
differences, compared with the case of the scalar field, are, apart from the appearance of 
the projection operator, the absence of the factor (l/2g + ) in the integrand of (13) and 
the negative sign of the second term in the fermionic case. They, however, compensate, 
and the standard manipulations to factor out the exponential give rise to the factor 
[0(q + ) + 9(— q + )] which may be interpreted as unity in the distribution theory sense, 
parallel to what we find in the derivation of the scalar field propagator on the LF. The 
straightforward use of the integral representation 2iri 9{r) e~ tpT = J dX e^~ lXr ' /(p — A — 
ie) of 0(±t), together with the standard manipulations in the second term to factor 
out the exponential, leads to 

< 0\T(^ + (x)^ + (0))\0 >= -— / d*q V \ e-*-', (14) 

(27t) 4 J (q z — m z + ie) 

where for convenience we have renamed the dummy integration variable originating 
from the integral representation of 9(±t) as q~ and d A q = d 2 q ± dq + dq~ with all integra- 
tions ranging from — oo to oo. It may also be derived straightforwardly by functional in- 
tegral; we do have to take care of the second class constraint in the measure. We have re- 
stored the color index as well. The propagator (14) contains no instantaneous term M, 
and the integrand factor may also be expressed as ~ [A + (jq + m)A~ /(q 2 — m 2 + ie)] 7 . 
We verify that the propagator satisfies the equation for the Green's function corre- 
sponding to the equation of motion of ip + , (8). 



4 Gauge field propagator 

We now turn our attention to the gauge field propagator on the LF in the covariant 
gauges. The relevant quadratic terms in the Lagrangian density are 

i [F\_F\_ + 2F a +± F a _ ± - F\ 2 F\ 2 ] + +B a 9,i ffi " + ^B a B a + td^xi'd^ (15) 



It will be sufficient to study the free abelian gauge theory described by the following 
Lagrangian density 

1 f , „ xO , „ s O _ „ „ 1 ,, , . , . , .Is £ „2 



2 



(F + _) 2 - (F 12 ) 2 + 2F +± F_ ± \ + B(d + A_ + d_A + + d ± A x ) + ^B\ (16) 



where for the abelian theory F^ = (9 M A,— <9j,A M ). The canonical momenta are ir + = 0, 
n b = 0, n 1 - = F_j_, tx~ = F + _ + B and the canonical Hamiltonian density is found to 
be 

H c =^(7t-) 2 + ^(F 12 ) 2 -A + (d_ir- + d ± 7r ± -2d_B)-B(n- + d ± A ± ) + ±(l-0B 2 (17) 



In the Dirac procedure the primary constraints [I I] are n + w 0, 7r^ w and r\ = 
it 1 - — d_Aj_ + d±A_ m 0, where _L= 1, 2 and ~ stands for weak equality relation ]14| . 
We now require the persistency in r of these constraints employing the preliminary 
Hamiltonian, which is obtained by adding to the canonical Hamiltonian the primary 
constraints multiplied by the Lagrange multiplier fields. We assume the standard 
Poisson brackets for the dynamical variables in the computation for obtaining the 
Hamilton's equations of motion. We are led to the following two secondary constraints 

$ = (9_tt- + diTT 1 - 2(9_5 w 0, 
^f = n- +2d^A + + d ± A ± -(1-£)B w 0. (18) 

The Hamiltonian is next enlarged by including these additional constraints as well. 
The procedure is repeated. No more constraints are seen to arise, and we are left only 
with the equations which would determine the Lagrange multiplier fields. According 
to the Dirac procedure Jbl , we now go over from the standard Poisson brackets to the 



modified Poisson brackets, called frequently Dirac brackets, such that inside them we 
are able to substitute the above constraints as strong equality relations (e.g., by = in 
place of ~). The equal-r Dirac bracket {f(x),g(y)}o which carries this property is 
constructed straightforwardly. Hamilton's equations now employ the Dirac brackets 
and the phase space constraints n + = 0, ttb = 0, rj = 0, $ = 0, and \I/ = then 
effectively reduce the (extended) Hamiltonian. In the covariant Feynman gauge with 
£ = 1 the free Hamiltonian takes the simple form 



Hi 



LF 



- f d 2 x ± dx~ g^A,, d ± d ± A v . (19) 



The theory is canonically quantized through the correspondence i{f(x),g(y)}n — > 
[f(x),g(y)], the commutator among the corresponding operators. The equal-r com- 
mutators of the gauge field are found to be [A^(x), A u (y)] x+ +=t = —ig IMU K(x,y) 

8 



where K(x,y) = — (l/4)e(x~ — y~)5(x ± — y 1 -) is nonlocal in the longitudinal coordi- 
nate. The transverse components of the gauge field have the physical LF commutators 
|y4_i_(x), A±r(y)] T = i8j_ t ±i K(x, y) , while for the ± components we have only the mixed 
commutator nonvanishing [A + (x), A_(y)] T = —iK(x,y), it has a negative sign which 
indicates the presence of unphysical degrees of freedom in Feynman gauge. For £ 7^ 1 
the commutator, for example, of A± with A± is found to be nonvanishing. 

The Heisenberg equations of motion lead to OA^ = for all the components, and 
consequently the Fourier transform of the free gauge field over the complete set of plane 
wave solutions takes the following form on the LF 



A»(x) = , l f d 2 k L dk + ^J=e^ x \k) \a {X] (k + ,k L )e~ ik - x + a) x Ak + ,k L )e ik 

(20) 
where e^ x \k), A = —,+,1,2 label the set of four linearly independent polarization 
four- vectors. In the front form theory the two transverse (physical) polarization vector 
are space-like while the others are null four- vectors. For a fixed k^ = (k°, k), where A; = 
I A; |, we may construct them as follows: e^ = (1, k/k )/\/2 , e^~> = (l,—k/k°)/\^2, 
e^ 1 ' = (0, e(k;l)) , and e^ = (0, e*(/c; 2)). Here (0,1,2,3) components are specified 
for convenience while e(k;l), e(k;2) and k/\k\ constitute the usual orthonormal set 
of 3- vectors. The polarization vectors are orthonormal: g^ v e^^{k)e^ u {k) = g Xa 
and satisfy the completeness relation: <7.\<7 e >(&) e i>(&) = g^ v . The field commu- 
tation relations for the gauge field found above are shown to be satisfied if we assume, 
oarallel to the discussion in the fermionic case, the canonical commutation relations: 
'a {x) (k + ,k ± ),a\ a) (k' + ,k' ± )] = -g Xa 5(k + - k' + ) S 2 (k ± - k' x ). We note that the op- 
erators a( ) = (a,( + ) + 0(_))/v2 and a^ = (a(+) — fl(_))/v / 2 obey the usual canonical 
commutation relations except that in the case of O(o) a negative sign is obtained. The 
discussion of the Gupta-Bleuler consistency condition then becomes parallel to that 
in the conventional equal-time treatment |y| of the theory. The Feynman gauge free 
gauge field propagator on the LF can be derived straightforwardly using (20) 

< 0\T(A%(x)A\(0))\0 >= ^ J d'k e~ ik -* ^- (21) 

where we have restored the gauge index a and used a(\) a ( y k + ,k ± ),a\ a )b( y k' + ,k' ± ) 
= _ 9\ababb(k + — k' + )5 2 ( y k ± — k'-*-). The free ghost-antighost propagators for the an- 
ticommuting ghost fields may likewise be derived; they are also causal and given by 
(V '(2tt) 4 )<W (k 2 + ie).) 



5 Illustrations 

The Hamiltonian density in Feynman gauge is 

n_/LF nj i nj 

rt = n.Q + ttint 

$\ [^'(n 1 ^ " m) + g^A^j] 4>J + 9 -f abc {d^A% - d u A%)A b »A c » 
+ i-f^f^A\A\A^A dv + d»(x a Kx a + gF^rh'A^ (22) 

where if) J is given in (4), we have set a/2% — (xi + %), V%X — (xi — %), and in 
(22) the cubic and higher order terms belong to 7i ini which is also understood to be 
normal ordered. It is worth remarking that despite the presence of the longitudinal 
operators a± and a^± in the fields Ap, there are no non-zero matrix elements involving 
these quanta as external lines in view of the commmutation relations of these operators 
as discussed in the previous section. 

The perturbation theory expansion in the interaction representation where we time 
order with respect to the LF time r is built following the Dyson- Wick [|U|] procedure. 
We will illustrate it in our context through some explicit computations, for simplicity, 
in QED where U(x\A_) = exp{—ief x du~ A^.(t,u~ ,x )} and D^ = (d^ — ieA^). We 
observe from (4) and (22) that a seagull term of the order e 2 is present in the interaction 
Hamiltonian at the tree level as in case of the scalar field QED. 

Consider Electron-Muon scattering. The contribution coming from the longitudinal 
components of the gauge field, viz, A + and A_, which carry a nonvanishing mixed 
propagator, is contained in the following normal-ordered product of the Wick decom- 
position of the second order term in the perturbation theory expansion of the S-matrix 

e 2 J d 4 x 1 d A x 2 : ^ e+ (x 1 )if) e+ (x 1 )A+(x 1 ) i)\ + (x 2 )(m + % r ) 

- x / dy 2 'e(x 2 ~ - y 2 ~){ / du 2 ~ A„(u 2 )}(m - i fi T )ipp + (y 2 ) : (23) 

where overdots indicate |0| the pair of fields contracted, x 2 = (r 2 ,x^",x^), u 2 = 
(r 2 ,wj,x^), y 2 = (r 2 ,yJ,X2"), X\ = {ri,Xi,xj;), and p T = j ± d±. Here we write 
J e(x — y)dy/2 in place of d x ~ and suppress the symmetry factors for convenience. 
Making use of the Fourier transforms of the fields and the gauge field propagator the 
matrix element is written down by simple inspection 



e 2 



(m 


+ />'/) 


(m 


-A T ) 




~¥ + 2 




-ip+ 2 



u\+{p'l)Ue+{Pl) U\+(p 2 ) —T — U^{p : 



-Tj- (24) 
q 2 + xe 



10 



where g M = {p'^ 2 — p M i). It may be readily rewritten in view of (10) and the simplifi- 
cations following from the mass shell conditions on the external lines to give 



e 2 



U e (Plh +U e(Pl) uM)l U »(P2)] o , + ■ ( 25 ) 



q 2 + it 



The graph with the A + and A_ interchanged gives rise to a similar expression with 

g + _ — > g h while 7 ± — > 7 T . There are four contributions arising from the virtual 

propagation of the transverse components of the gauge field. They add up to 



e 2 



«e(pi)7 ± «e(Pl) U^p'^'u^)] f 1 ' (26) 

J q z + te 
resulting in the following complete matrix element to the second order 

e 2 [« e (ft)7\(Pi)«,(p' 2 )7X(p 2 )] -=^-- (27) 

which agrees with the one obtained in the conventional equal-time formulation. 

Consider next the computation of Electron Self-Energy. The contribution from the 
longitudinal components arises from 



/oo ] 
-dy 2 ~e{x 2 
-00 2i 



2/2 



x-z" 



T\j\ (.. w,t 



{ / du<f A_(u 2 )}(m - i fa 1 )^ + (y 2 )^+(x 1 )il; + (x 1 )A + (xi) : (2 



in 



leading to 



2 f , A uW(p)[>y-(m+ A T )1 + W 8) (P) ( „ x f9Q x 

J [{p — q) z + %e\ {q z — m z + te) 

The graph with the A + and A_ interchanged gives rise to a similar expression with 
<7+_ — > g_+ while 7 ± — ► 7 T . The matrix elements following from the four graphs 
corresponding to the exchange of the ( photon ) fields A± and A 2 is also written down by 
simple inspection. As in the earlier case the expressions get simplified in virtue of (10) 
and acquire the covariant form encountered in the conventional covariant perturbation 
theory. The complete matrix element is found to be 

J [{p — q) z + ie\{q z — m z + it) 

where g M = ((m 2 +q- L q- L )/(2q + ), q + , q x ) and the integration measure is d 4 q = d 2 q ± dq + dq~ 
Considering that the integrand has the pole at q 2 —m 2 RiOwe may regard the expression 
obtained on the LF to be effectively identical to the one obtained in the conventional 
covariant theory framework. The contribution coming from the seagull term at the 
tree level vanishes if the dimensional regularization is used. 
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The calculation of the tree graphs for Compton scattering, 7 + e — > 7 + e, is tedious 
but straightforward. The results on the LF are shown to be in agreement with the 
conventional covariant theory one. We remark that on the LF the tree level seagull 
term dominates the (classical) Thomson formula for the scattering at the vanishingly 
small photon energies similar to the case of the QED with the scalar fields. In the 
covariant gauges on the LF all the relevant field propagators in momentum space 
are causal. Employing (10) and (12) it is rather straightforward to rewrite the final 
result in manifestly covariant form. The computation to higher orders in the Dyson- 
Wick perturbation expansion can be carried out straightforwardly; the nonlocality 
of the interaction, arising from (4), does require some extra effort to handle, but it 
becomes easier to control the rotational invariance and compute the loop integrals in 
the traditional fashion. 



6 Conclusions 

The LF Dyson- Wick perturbation theory expansion based on the LF-time-ordering has 
a number of advantages for computing the high order corrections in the front form QCD 
theory. The covariant Feynman (or Landau) gauge may be adopted with advantage on 
the LF where all the relevant propagators take (except for the numerical projector in 
the dynamical fermion propagator) the covariant causal forms permitting us to employ 



the usual power counting rules. The Hamiltonian version |26| can be implemented 



in DLCQ |T9|. In the approximation of no retardation, the Coulomb interaction is 



recovered. The propagator of the independent component of the fermion field on the 
LF has been shown to be simpler than in the conventional theory. The momentum 



space integrals in the front form theory may be converted |2T] to the Euclidean space 
integrals which then permit us to employ, for example, dimensional regularization. The 
illustrations given here demonstrate the agreement of the LF quantized theory results 
with the conventional covariant theory ones. The expression (10) for the LF spinor 
proved to be quite useful here. The fact that in the front form theory the classical 
Thomson scattering limit is obtained from a seagull term at the tree level is significant 
since, it seems difficult to build on the LF a systematic procedure to obtain semiclassical 
approximation. 

It is worth remarking also that we have made an ad hoc choice of only one (of the 
family) of the characteristic LF hyperplanes, x + = const., in order to quantize the 
theory. The conclusions here confirm the conjecture 0] made earlier on the irrelevance 
in the quantized theory of the fact that the hyperplanes x ± = constitute characteristic 
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surfaces of hyperbolic partial differential equation. The discussion given in this paper 
may clearly be repeated for the case of the non-covariant gauges such as <9_A_ = 0. 
It will appreciably simplify the expression of ip- given in (4) , but now in view of the 
form of the gauge field propagator, we will require the well known prescriptions [27| in 
order to recover on the LF the conventional covariant theory results. 
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